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Abstract. We present results of our 1.4 GHz image stacking analysis of mass- 
selected galaxies in the COSMOS field. From the resulting median radio continuum 
flux density we have determined the evolution of the average star formation rate (SFR) 
of galaxies as a function of stellar mass, unbiased from effects of dust but also source 
confusion due to the 1 .5" angular resolution achieved by the VLA. We find a power-law 
relation between specific SFR (SSFR) and stellar mass for star forming galaxies out to 
Z = 3. While higher mass systems exhibit lower SSFRs at any epoch, no differential, 
more rapid evolution of high mass galaxies is evident at least out to z ~ 1.5 where 
our conclusions are most robust. Utilizing measured mass functions of star forming 
systems, the characteristic stellar mass for galaxies contributing most to the comoving 
SFR density appears not to evolve. These findings hence challenge 'downsizing' sce- 
narios in which star formation has gradually shifted towards lower mass systems with 
cosmic time. Our analysis constitutes the to-date best determination of the cosmic star 
formation history (CSFH) since z = 3 and yields indirect evidence for a rapid decline 
of the global mass density of molecular gas with time. 



1. Introduction 

A complete census of the global cosmic star formation activity at a given cosmic epoch 
is usually obtained by measuring/ extrapolating the total star formation rate (SFR) per 
unit comoving volume. The time evolution of this SFR density (SFRD) is also known 
as the cosmic star formation history (CSFH). Over the last years several measurements 
revealed a steep decline of the SFRD sinc e a redshift of z ~ 1 and like ly a shallower evo- 
lutionary behavior at higher redshifts (see lHopkins & Beacomll2006l for a compilation). 
Despite the to-date comparatively little observational evidence into the SFRD evolution 
at high z - and in particular into the fraction of star formation obscured by dust - it is a 
widely accepted scenario that the CSFH shows a maximum around z = 2 - 3. 

A number of recent studies revealed a surprisingly tight correlation b etween the to 



tal SFR of normal star forming galaxies and their stellar mass content (e.g. lNoeske et al. 



20071 : lElbaz et al.l l2007. for measurements out to z ~ 1). The localized process of star 



formation is consequently linked to a global galaxy property, suggesting that star for- 
mation histories of individual galaxies are Ukely n ot dominated by stochastic processes 



such as galaxy mergers (e.g. iNoeske et al.ll2007b . The exact evolutionary properties 



of this relation are subject to ongoing debates and insights based on a variety of SFR 
tracers are needed. 
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Clearly, multi-wavelength look-back surveys are key to reveal how galaxies build 
up their stellar mass and how their SFRs evolve over a wide range of reds hift. The 



panch romatic datasets from the 2 deg cosmic evolution survey (COSMOS: see lScoville et al. 



2OO7L for an overview), in particular, have provided very accurate photometri c redshifts 



(photo-zs) and stellar masses for a large mass-selected sample of galaxies (Ilber t et alj 
I201Q) . In this work we use this Spitzer/IRAC (mAB(3.6 /im) < 23.9) selected sample 
of ~ 114,000 galaxies in combination with high angular resolution Very Large Array 
(VLA) radio continuum data to measure the evolution of the average SFR as a function 
of stellar mass out to z ~ 3. SFRs derived from radio continuum emission require only 
a simple K-correction and no correction for du st extinction. AH m ethods, results and 
Figures presented here are discussed in detail in iKarim et al.l (l201lh . 



2. The evolution of the specific star formation rate as a function of stellar mass 



Stellar masses have been estimate d for all sources via fits of stellar population synthesis 
models (iBruzual & Chariot 2003, BC03) to the observed spectral energy distributions 
(SEDs), assuming a iChabrier (2003.) init i al ma ss function, exponentially declining star 
formation histories and a ICalzetti et al.l (l2000h dust exti nction law (detail s on the es- 
timation of stellar masses and photo- z's are found in lllbertetal.l[2010l) . From the 
restframe intrinsic - i.e. dust unextincted - (NUV - r"'^)temp < 3.5 colors, obtained from 
the BC03 SED fits, we have separated star forming sources from quiescent ones within 
our sample. 

We have obtained average 1.4 GHz radio flux densities separately for both all and 
only star forming galaxies fr om median stacks of cu t out im ages of the VLA-COSMOS 



map (resolution 1.5" x 1.4": ISchinnerer et al.ll2007l l2010l) centered at the positions of 



the 3.6 selected sources, after binning each galaxy sample in redshift and mass. 
We also exc luded radio-AGN c andidates from the sample - applying various criteria 
discussed in iKarim et al.l (120 111) - prior to stacking. The resulting total flux densities 
have been converted into SFRs using the calibration of the radio-infrared correlation in 
the local universe (Bell 2003). 

The average stellar mass normalized SFR (specific SFR/SSFR) is best described 
by separable power-laws in mass and redshift for both all (massive) galaxies (Fig. [T] 
left panels) and star forming ones (Fig. [T] right panels). At z > 1-5 the average SSFR 
of low mass star forming galaxies tends to deviate from the mass-uniform evolution- 
ary trend seen at lower redshifts. There is thus tentative evidence for an upper limit to 
the average SSFR at just about the inverse of a dyn amical timescale of typical local but 
also higher z disk galaxies (e.g. lDaddi et al.ll2010al) . While their dynamical and free-fall 
times are approximately equal, normal galaxies are assumed to accrete thei r gas from 



their su rroundings and to process it into stars at a constant efficiency (e.g. Daddi et al. 



2010bh . The upper SSFR-limit we propose might hence represent an effective gas ac- 
cretion threshold. 



3. The characteristic mass of star formation 



While star forming sources exhibit this remarkably shape-invariant SSFR-mass relation 
throughout cosmic time in our study, it is worth noting that their stel lar mass function s 
are found to show a similar constancy of their Schechter profile (e.g. lllbert et al.ll20Tol) . 
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Figure 1. Radio stacking based measurement of the SSFR as a function of stel- 
lar mass at 0.2 < z < 3.0 (upper panels) and redshift evolution of the SSFR in 
logarithmic stellar mass bins (lower panels) for our entire galaxy sample (left) and 
star forming systems only (right). Full symbols depict samples that are regarded as 
representative for the underlying galaxy population, open symbols represent upper 
SSFR limits. For representative data points both the mass dependence of the av- 
erage SSFR at a given epoch as well as its redshift evolution at fixed mass follow 
power-laws (dashed lines). Quiescent galaxies (Included in the left panels) are pref- 
erentially massive and less frequent at higher z as seen in (i) deviations from the 
high-mass end power-law which itself (ii) becomes gradually slightly shallower with 
redshift. For star forming systems only (right panels) the power-law extends over 
the entire mass range and all mass-bins evolve uniformly. A separable power-law 
SSFR(M,, z) oc f(M») X g(z) - M;° °'* x (1 -H z)^-^ hence is a good description for the 
average SSFR of star forming galaxies at least out to z ~ 1.5. The horizontal band 
(upper panels) sketches the inverse dynamical timescale, potentially representing an 
upper bound to the average SSFR (see text) and therefore leading to deviations also 
for high-z star forming sources. The black long-dashed line gives the mass-doubling 
limit above which galaxies are able to double their mass until z = assuming a con- 
stant SFR. The black dashed-dotted line depicts the inverse age of the universe at any 
given redshift and hence makes measured SSFRs co mparable to the past av erage star 
formation activity. The SED-derived measurement of lMagdis et alJ (1201 Ol) for LBGs 
at z ~ 3 with log(M,) ~ 10 is shown as a filled star 
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Figure 2. The distribution of the SFRD with respect to stellar mass at different 
epochs out to z ~ 3. The data points in each panel are the product of the observed 
number densities of star forming galaxies with the average (stacking- based) radio 
SFRs. Leftwards of the dashed vertical lines - in the mass range where our data 
are not fully representative for the undelying galaxy populat ion - number den sities 
have been corrected using the corresponding mass functions dllbert et al.ll2010l) . Un- 
corrected data are shown as open circles. These corrections are generally small and 
needed only at z > 1.5. In each panel we overplot the Schechter function that re- 
sults from multiplying the best-fit radio derived SFR-mass relation at a given epoch 
with the corresponding mass function fo r SF galaxies. Literature results based on 
the [O IIU3727 Une (Gil banket al.ll201d) as well as UV emission dCowie & Barged 
12008) agree well with the trends in our data. Dashed blue lines show the distribution 
obtained if an upper limit to the average SSFR at lower masses is assumed. 



As both phenomena involve a power-law in stellar mass, the product of the S SFR- 
mass relation and the mass function at a given cosmic epoch also yields a Schechter 
function for the mass distribution of the comoving SFR density (SFRD). In Fig. [2] 
we plot the resulting Schechter SFRD functions along with our direct, radio stacking 
based, measurements of the SFRD for each bin in mass and redshift and with literature 
results from alternative SFR diagnostics at different wavelengths. While the potential 
upper SSFR-limit discussed above might alter the low-mass end shapes significantly 
only at the highest redshifts, the galaxy mass at which the function peaks at any epoch 
is Mh. ~ lO'"-^ Mq, regardless of the assumption of an SSFR-limit. The existence of 
this characteristic mass of star forming galaxies challenges a 'downsizing' paradigm in 
which the dominant contributors to the CSFH have been massive sources in the cosmic 
past and low-mass galaxies at present times. It does, however, not rule out that low- 
mass sources gain more relative importance over time. 
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Figure 3. The dust-unbiased CSFH out to z = 3. Black circles represent the sum 
of data points the corresponding redshift-bin panel of Fig. |2] and hence a direct 
- stacking-based - measurement of the SFRD down to the limiting mass at each 
epoch. Red filled circles correspond to the 'total' SFRD at each epoch, obtained by 
integrating the Schechter-function fit (Fig. |2]i down to M, = 10^ Mq and assuming 
an upper limit to the average SSFR. The redshift evolution can be described by a 
broken power-law (dashed lines) that results from the joint (non-)evolution of the 
star forming stellar mass density and the evolution of the SFR-mass relation. Re- 
cent mid- to far- IR measurements of SFRDs bet ween < z < 2.5 (Spitzer/MIPS: 
iRodighiero et al.ii2010: Herschel/R4CS: GruDDioni et al.'2010l) are shown along with 
our data and the 3(T envelope from the Hopkins & Beacom ( 2006) compilation. Open 
symbols denote lower limits. Note the remarkable agreement of the far-IR- and 
radio-based data at all z- 



4. The dust unbiased cosmic star formation history 



The sum of our direct radio stacking based measurements of the SFRD for individual 
mass bins yields lower limits to the total SFRD at any redshift. At a given redshift the 
Schechter SFRD function enables us to obtain the remainder of the total SFRD, inte- 
grated below our observational mass limits. For this integration we assume the upper 
bound to the average SSFR motivated above. The resulting radio-based - i.e. dust- 
extinction unaffected - CSFH out to z ~ 3 is shown in Fig. [3] along with independent 
literature measurements. Since z ~ I the stella r mass density of star forming galaxies is 
approximately constant (e.g. Ilb ert et al1l2010h and also the relative fractions of objects 
at different masses within this population are conserved due to the shape-invariance of 
their mass function (see above). In this redshift range the CSFH is therefore entirely 
controlled by the mass-uniform evolution of the average SFR. Given a constant star 
formation efficiency, a strong and global decline in the mass density of molecular gas 
thus entirely explains the observed decrease of the integrated SFRD with cosmic time. 
Towards earlier epochs it is predominantly the - again mass- uniform - decrease in star 
forming stellar mass density that leads to a shallower evolutionary behavior. The CSFH 
since z ~ 3 is hence well described by a broken power-law (Fig. |3]l. 
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5. Summary and conclusions 

Using a median image stacking technique of 1.4 GHz radio continuum emission and an 
unprecedentedly rich sample of galaxies selected at 3.6 yum with panchromatic (FUV to 
mid-IR) ancillary data in the COSMOS field we have measured stellar mass-dependent 
average (specific) stai" formation rates ((S)SFRs) in the redshift range 0.2 < z < 3. We 
have applied various criteria to minimize contaminating radio flux from potential active 
galactic nuclei and have used the template-based, intrinsic rest-frame (NUV - r''')temp 
color from SED-fits in the NUV-mid-IR in order to sep arate star forming s ources from 



quiescent systems. Our findings - presented in detail in lKarim et al.l ( 20Tll ) - are: 



• The evolution of the average (S)SFR of stellar mass selected galaxies is approx- 
imately mass-independent at least out to z ~ 1.5 where the depth of our data 
allows the most robust conclusions. A power-law relation between (S)SFR and 
mass represents well the data at the high mass end. 

• The (S)SFR of star forming sources at z < 3 also evolves basically independent of 
mass while a power-law dependence of (S)SFR and mass generally extends over 
all masses probed. Higher mass galaxies hence have lower SSFRs, regardless if 
quiescent galaxies are included in the analysis or not. 

• There is tentative evidence for an upper bound to the average SSFR preventing a 
further growth of the SSFR with z, starting at the low-mass end at z ~ 1.5-2. 

• The stellar mass distribution function of the comoving SFR density at any epoch 
is well described by a Schechter function that reveals a constant characteristic 
mass of star formation since z ~ 3. Its integral yields the total SFRD at any z. 
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